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ELF/VLF  Electromagnetic  Detection  and  Characterization 
of  Deeply  Buried  Targets 


This  is  the  Final  Report  for  the  SBIR  Phase  II  effort  entitled  "ELF/VLF  Electromagnetic 
Detection  and  Characterization  of  Deeply  Buried  Targets"  sponsored  by  the  Defense 
Advanced  Research  Projects  Agency.  The  various  electromagnetic  calculational  tools 
coded  as  MATLAB  routines  delivered  with  the  report  constitute  the  goal  of  our  effort. 
The  report  is  in  five  parts:  (1)  a  description  of  fhe  Sommerfeld  roufines;  (2)  a  description 
of  the  menu  interface  routines;  (3)  a  description  of  fhe  Hill  /  King  buried  wire  routines; 
(4)  an  appendix  of  the  Sommerfeld  infegrals  used;  (5)  an  appendix  describing  fhe  CD- 
ROM  containing  the  MATLAB  codes. 

Sommerfeld  Calculation  Routines  for  a  Two-Layer  Earth 

The  configuration  for  fhe  Sommerfeld  roufines  is  sketched  below: 


z 


To  calculate  the  electric  or  magnetic  fields  due  to  a  dipole  at  a  set  of  points, 
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MATLAB  function  routines  are  used:  getElecDipoleField  or  getMagDipoleField; 
these  functions  are  called  in  MATLAB  as: 


[EField , HField]  =  ... 

getElecDipoleFieldCEFieldOn , HFieldOn , ILdipMom , . . . 

RSource Point, RE vat PointSjfreqVal, . . . 

zLayerl , eDielecl , Conductl , zLayerZ ,eDielec2,Conduct2); 


[EField , HField]  =  ... 

getMagDipoleFieldCEFieldOn , HFieldOn , lAdipMom , . . . 

RSource Point, RE valPointSjfreqVal,  .  .  . 

zLayerl , eDielecl , Conductl , zLayer2 ,eDielec2,Conduct2); 

The  description  for  most  of  the  arguments  is  the  same  for  both  routines;  the  routines' 
input  arguments  differ  only  in  the  definition  of  the  dipole  moments. 


Below  are  the  input 
EFieldOn 

HFieldOn 

RSourcePoint 


RFvalPoints 


freqVal 


argument  descriptions: 

-  Logical  variable:  ON  (=1)  ->  Calculate  the  E  Field; 

OFF  (=0)  ->  Don't  Calculate  the  E  Field 

-  Logical  variable:  ON  (=1)  ->  Calculate  the  H  Field; 

OFF  (=0)  ->  Don't  Calculate  the  H  Field 

-  Location  of  the  dipole  source  in  meters  for  user  defined 
Cartesian  coordinates;  RSourcePoint(l,3), 

i.e.  RSourcePoint  =  [xsrc,ysrc,zsrc] 

-  Points  to  evaluate  the  E  and/ or  H  Fields  in  meters;  uses  the  same 
(user  defined)  Cartesian  coordinate  system  as  RSourcePoint 
(above);  REvalPoints(l:N,3),  where  N  is  the  number  of  eval  points 
i.e.  REvalPoints(l,:)  =  [xEv(l),yEv(l),zEv(l)], 

REvalPoints(2,:)  =  [xEv(2),yEv(2),zEv(2)], 

REvalPoints(3,:)  =  [xEv(3),yEv(3),zEv(3)], 

. etc. 

N.B.  For  a  single  call  of  either  routine  the  REvalPoints  selected 
must  all  be  in  a  single  medium;  if  the  E/H  Field  values  are  desired 
in  both  Air  and  Region  1,  this  may  be  accomplished  in  two 
sequential  calls.  Also  note:  Fields  cannot  be  calculated  in  Region  2, 
as  this  is  intended  to  model  a  substrate  below  the  air  and  ground 
regions  of  sources  and  observation  points 

-  Frequency  in  Hz 
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zLayerl 

eDielecl 

Conductl 

zLayer2 

eDielec2 

Conduct2 

ILdipMom 

lAdipMom 


-  z  value  [meters]  of  Air  <— >  Region  1  interface  in  fhe  user  defined 
Carfesian  coordinate  system 

-  Relative  dielectric  constant  of  Region  1  medium 

-  Conductivity  of  Region  1  medium  [mhos  /  m] 

-  z  value  [mefers]  of  Region  1  <— >  Region  2  interface  in  fhe  user 
defined  Cartesian  coordinate  system 

-  Relative  dielectric  constant  of  Region  2  medium 

-  Conductivity  of  Region  2  medium  [mhos  /  m] 

-  Electric  dipole  moment  [amp-mj;  used  only  for  fhe 
function  getElecDipoleField 

-  Magnetic  dipole  moment  [amp-m^j;  used  only  for  fhe 
funcfion  getMagDipoleEield 


Below  are  fhe 
EEield 


output  value  descriptions: 

-  The  value  of  fhe  E  Field  due  to  the  dipole  at  the  specified 
evaluation  points  (REvalPoints)  in  volt/ m.  The  E  Field  is  given  in 
Cartesian  coordinates.  So, 

EField(l,:)  =  [Ex,Ey,Ez]  at  REvalPoints(l,:) 

EField(2,:)  =  [Ex,Ey,Ez]  at  REvalPoints(2,:) 

EField(3,:)  =  [Ex,Ey,Ez]  at  REvalPoints(3,:) 

. etc. 


HField  -  The  value  of  fhe  H  Field  due  to  the  dipole  at  the  specified 

evaluation  points  (REvalPoints)  in  amp  /  m.  The  H  Field  is  given  in 
Cartesian  coordinates.  The  array  structure  and  correspondence 
between  HField  and  REvalPoints  is  the  same  as  for  EField  above. 


There  are  a  number  of  MATEAB  functions  required  by  getElecDipoleField  and 
getMagDipoleField;  fhese  functions  are  supplied  on  the  CD-ROM  and  they  must  be 
kept  in  the  same  folder  as  getElecDipoleField  and  getMagDipoleField. 

Example:  Simple  Buried  HED  Calculation 

As  a  first  example  we  give  below  the  code  for  calculation  of  an  HED  in  Region  1  at  a 
depth  of  0.5  meter  below  the  surface  at  2  kHz.  The  relative  dielectric  constant  of  Region 
1  is  10  and  fhe  conductivity  is  0.01  mhos/ m;  the  relative  dielectric  constant  of  Region  2 
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is  20  and  the  conductivity  is  0.05  mhos/ m.  The  interface  between  the  Air  and  Region  1 
corresponds  to  z  =  0  m;  the  interface  befween  fhe  Region  1  and  Region  2  corresponds 
fo  z  =  -100  m.  We  assume  fhe  elecfric  dipole  sfrengfh  is  0.5  amp-m.  The  dipole  is 
cenfered  af  x  =  0  m,  y  =  0  m,  and  z  =  -0.5  m  in  our  coordinafe  sysfem  and  fhe  dipole 
vecfor  is  along  fhe  y-axis.  The  elecfric  and  magnetic  fields  are  evaluafed  on  a  square 
grid  30  m  X  30  m  a  disfance  of  1.0  m  above  fhe  surface  af  2  mefer  infervals.  The  resulfs 
are  ploffed  for  each  field  componenf  using  fhe  MATLAB  pcolor  command. 

The  configuration  for  fhis  example  is  skefched  below. 


02=0.05  mhos/m 


Region  2 


The  MATLAB  code  for  fhis  calculation  is: 


%  Set  flags  to  calculate  both  E  Field  and  H  Field 

EFieldOn  =  l;HField0n  =  1; 

%  Frequency  =  2  kFIz  and  set  dielectric/conductivity  for  Region  1  &  2 

freqVal  =  2000; 

ELF/VLF  SBIR  Phase  II  Final  Report  CODAR  Ocean  Sensors,  Ftd.  Page  4  of  42 


eDielecl  =  10;Conductl  =  0.01; 
eDielecZ  =  20;Conduct2  =  0.05; 

%  Region  1  <-->  Air  at  z  =  0  m;  Region  2  <-->  Region  1  at  z  =  -100  m; 

zLayerl  =  0;zLayer2  =  -100; 

%  Electric  dipole  moment  =0.5  Amp-m,  along  y-axis 

IL  =  0.5; 

dipMom  =  IL. *[0,1,0]; 

%  Electric  dipole  moment  at  x  =  0  m,  y  =  0  m,  z  =  -0.5  m 

%  Since  the  Region  1  <-->  Air  interface  is  at  z  =  0  m  (zLayerl  =  0  m), 

%  this  places  the  dipole  o.5  m  below  the  earth 

RSourcePoint  =  [0,0, -0.5]; 

%  Set  up  a  square  grid  of  30  m  x  30  m  at  a  distance  of  1.0  m  above  the 

%  Region  1  <-->  Air  interface  at  2  meter  intervals  on  which  to 

%  evaluate  the  electric  and  magnetic  fields 

XX  =  CC-30):2;30). ' ;yy  =  CC-30):2;30). ' ;zz  =  1.0; 

xxTemp  =  xx*ones(sizeCyy . ' )) ;yyTemp  =  onesCsizeCxx))*Cyy . ' ) ; 

xxEval  =  xxTempC ; ) ;yyEval  =  yyTemp(;); 

npts  =  max(sizeCxxEval)) ; 

REvalPoints  =  zerosCnpts , 3) ; 

REvalPointsC : , 1)  =  xxEval ; REvalPointsC : , 2)  =  yyEval; 

REvalPointsC : , 3)  =  zz . *onesCnpts , 1) ; 

%  Do  calculation  to  get  E/H  Fields  at  the  points  on  the  evaluation  grid 

[EField , HField]  =  getElecDipoleFieldCEFieldOn , FIFieldOn , dipMom ,  .  .  . 

RSource Point, RE valPoints,freqVal, . . . 

zLayerl , eDielecl , Conductl , z Layer 2 ,eDielec2,Conduct2); 


{  MATLAB  code  to 


do  plots  } 


Note  that  it  is  not  necessary  to  just  use  points  calculated  in  the  x-y  plane;  any  set  of 
points  is  acceptable,  as  long  as  all  the  points  in  a  single  call  are  located  in  a  single  layer  - 
either  Air  or  Region  1.  We  selected  the  set  of  poinfs  fo  be  in  fhe  x-y  plane  for  fhis 
example,  jusf  fo  show  one  useful  configuration. 

The  resulfs  of  fhis  calculation  are  shown  on  fhe  following  page  in  fhe  four  color 
infensify  plofs  for  fhe  fhree  E-Field  componenfs  and  fhe  magnitude  of  fhe  E-Field  along 
wifh  similar  plofs  for  fhe  H-Field.  We  have  supplied  fhe  complefe  MATLAB  code, 
ExampleHED.m,  including  fhe  plotting  sequences,  for  fhis  routine  on  fhe  CD-ROM. 
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Example:  Generic  Geophex  GEM  Sensor  and  a  Buried  Iron  Sphere 

For  the  second  example  we  examine  a  calculation  for  a  sensor  system  similar  to  the 
GEM  sensors  developed  by  Geophex,  Ltd.  of  Raleigh,  Norfh  Carolina  (our 
subconfracfor  in  fhis  efforf).  The  sensor  sysfem  and  measuremenf  configuration  are 
skefched  below: 


z 


For  ^  regions: 

More  defails  on  fhe  GEM  sensors  are  given  in  our  Phase  1  Einal  Reporf  and  are,  of 
course,  available  from  Geophex.  We  are  using  a  somewhaf  generic  version  of  fhe  GEM 
in  fhis  calculafion  and  fhe  resulfs  will  fherefore  be  represenfafive  in  a  general  sense,  buf 
are  nof  meanf  fo  be  a  precise  model  for  any  specific  version  of  fhe  GEM  sensor. 

In  fhis  example  a  0.5  m  radius  iron  sphere  is  buried  1.5  mefers  below  fhe  surface  and 
excifed  by  a  source  coil  af  one  end  of  fhe  sensor  sysfem.  This  induces  a  magnetic  dipole 
momenf  on  fhe  sphere,  which  reradiafes  and  is  sensed  by  fhe  receiver  coil  af  fhe  ofher 
end  of  fhe  sensor  sysfem.  The  receiver  also  picks  up  fhe  direcf  field  from  fhe  source  coil 
and  so  fhe  fofal  signal  is  a  sum  of  fhe  direcf  field  (fhe  field  if  fhere  sphere  were  nof 
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present)  and  the  perturbed  field  due  to  the  sphere.  Since  the  GEM  is  designed  to 
subtract  out  the  direct  signal,  the  effective  "measured"  signal  is  the  voltage  ratio  of  fhe 
perfurbed  signal  fo  fhe  direcf  signal;  fhis  is  usually  given  in  parfs-per-million  (PPM). 

The  calculation  is  done  af  1200  Hz;  fhe  relative  dielecfric  consfanf  of  Region  1  is  10  and 
fhe  conducfivify  is  0.01  mhos/ m;  fhe  relafive  dielecfric  consfanf  of  Region  2  is  50  and 
fhe  conducfivify  is  0.1  mhos/ m.  The  inferface  befween  fhe  Air  and  Region  1 
corresponds  fo  z  =  0  m;  fhe  inferface  befween  fhe  Region  1  and  Region  2  corresponds 
fo  z  =  -3  m.  We  assume  fhe  magnetic  dipole  sfrengfh  of  fhe  source  coil  is  3  amp-ml  The 
sphere  is  cenfered  af  x  =  0  m,  y  =  0  m,  and  z  =  -1.5  m  in  our  coordinafe  sysfem.  The  z 
componenf  of  fhe  perfurbed  and  direcf  H-Field  are  evaluafed  on  a  square  grid  6m  x  6m 
a  disfance  of  1.0  m  above  fhe  surface  af  ~0.2  mefer  infervals.  The  resulfs  are  plotted  for 
each  field  componenf  using  fhe  MATLAB  pcolor  command. 

The  MATLAB  routine  calls  a  function  getMagDipMomCond.  This  function  is  supplied  on 
fhe  CD-ROM.  The  function  calculafes  fhe  induced  magnefic  dipole  momenf  coefficienf 
for  a  conducting  sphere  for  a  given  permeabilify  and  dielecfric  consfanf.  The  relafion 
was  derived  by  Waif  {Geophysics,  16,  No.  4,  666-672,  October  1951  and  Radio  Science,  Vol.  3 
(New  Series),  No.  10, 1030-1034  October  1968). 

The  MATLAB  code  for  fhis  calculation  is: 

Somme rf el dGlobals ; 

%  Source  Loop  strength  3  amp*m^2  and  the  loop  is  held  to  be 

%  parallel  to  the  earth,  so  the  dipole  moment  is  along  the  z-axis 

IA_GEM  =  3; 

dipMomGEM  =  IA_GEM. * [0 ,0 , 1] ; 

%  Source  and  sensor  loops  are  separated  by  2  meters  along  the  y-axis 
LGEM  =  2; 

%  The  initial  position  of  the  source  is  1  meter  (z  =  1  meter)  above 

%  the  earth;  it  is  separated  by  LGEM  from  sensor  along  y-axis. 

RSourcePoint0  =  [0, -LGEM/2 , 1] ; 

REvalPoints0  =  [0,+LGEM/2,l] ; 

%  Set  up  measurement  grid. 

LL  =  3;NN  =  30; 

nx  =  NN;xMin  =  -LL;xMax  =  +LL;dx  =  2*LL/Cnx-1); 

XX  =  xMin:dx:xMax; 
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ny  =  NN;yMin  =  -LL;yMax  =  +LL;dy  =  2*LL/Cny-1); 
yy  =  (yMin ; dy :yMax) . ' ; 

freqVal  =  1200; 

wVal  =  2*pi*f reqVal ; 

eDielecl  =  10;Conductl  =  0.01; 
eDielec2  =  50;Conduct2  =  0.1; 

zLayerl  =  0;zLayer2  =  -3; 

%  Set  flags  to  calculate  only  the  H  Field 

EFieldOn  =  0;FIFieldOn  =  1; 

%  Get  Direct  Field  at  Sensor  for  each  of  the  measurement  position  points 
%  on  the  grid 
REvalPoints  =  zerosCl,3); 

FIFieldDirZ  =  zeros(ny ,  nx)  ; 
for  iy=l;ny 

for  ix=l;nx 

gridPoint  =  [xxCl , ix) ,yyCiy , 1) ,0] ; 

RSourcePoint  =  RSourcePoint0  +  gridPoint; 

REvalPoints  =  REvalPoints0  +  gridPoint; 
[EFieldDirTempjFIFieldDirTemp]  =  ... 

getMagDipoleFieldCEFieldOn  ,  FIFieldOn  ,  dipMomGEM ,  .  .  . 

RSource Point, RE valPoints,freqVal, . . . 
zLayerl , eDielecl , Conductl , z Layer 2 ,eDielec2,Conduct2); 
FIFieldDirZCiy ,  ix)  =  FIFieldDirTempCl,3)  ; 

end 

end 

%  Get  Field  at  Underground  Sphere  for  each  of  the  measurement  position 

%  points  on  the  grid.  Sphere  located  at  x=y=0,  z=-1.5  m 

REvaUG  =  [0,0, -1.5]; 

FIFieldUGX  =  zerosCny ,  nx)  ; 

HFieldUGY  =  HFieldUGX ; HFieldUGZ  =  HFieldUGX; 
for  iy=l;ny 

for  ix=l;nx 

gridPoint  =  [xxCl , ix) ,yyCiy , 1) ,0] ; 

RSourcePoint  =  RSourcePoint0  +  gridPoint; 

[EFieldUGTemp ,  FIFieldUGTemp]  =  ... 

getMagDipoleFieldCEFieldOn  ,  FIFieldOn  ,  dipMomGEM ,  .  .  . 

RSource Point, REvaUG, f reqVal , . . . 

zLayerl , eDielecl , Conductl , z Layer 2 ,eDielec2,Conduct2); 
HFieldUGXCiy,ix)  =  HFieldUGTemp(l , 1) ; 

HFieldUGYCiy,ix)  =  HFieldUGTempCl , 2) ; 

FIFieldUGZCiy,ix)  =  FIFieldUGTempCl , 3)  ; 

end 

end 

%  Get  Perturbed  Field  at  Sensor  for  each  of  the  measurement  position 

%  points  on  the  grid 

epsFreeSpace  =  8.85E-12;muFreeSpace  =  4E-7*pi; 

epsOut  =  epsFreeSpace*CeDielecl+i*Conductl/CwVal*epsFreeSpace)) ; 
muOut  =  muFreeSpace; 

sigiron  =  1.03E7;muIron  =  3*muFreeSpace ; 

epsSph  =  epsFreeSpace*(0.0+i*sigIron/CwVal*epsFreeSpace)) ; 
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muSph  =  mulron; 

ax  =  0.5; ay  =  0.5;az  =  0.5; 

[inducedDipMomCoef]  =  ... 

getMagDipMomCondCf reqVal ,epsOut,muOut,epsSph,muSph,ax,ay,az); 

RSourceUG  =  REvaUG; 

HFieldPertbZ  =  zerosCny , nx) ; 
for  iy=l;ny 

for  ix=l;nx 

gridPoint  =  [xxCl , ix) ,yyCiy , 1) ,0] ; 

REvalPoints  =  REvalPoints0  +  gridPoint; 

HFieldUG  =  [FIFieldUGXCiy , ix)  , FIFieldUGYCiy  ,ix),FIFieldUGZCiy,ix)]  ; 
inducedDipMom  =  inducedDipMomCoef .  *FIFieldUG ; 

[EFieldPertbTemp ,  FIFieldPertbTemp]  =  ... 

getMagDipoleFieldCEFieldOn  ,  FIFieldOn  ,  i  nducedDipMom  ,  .  .  . 
RSourceUG, RE val Points, freqVal, . . . 

zLayerl , eDielecl , Conductl , z Layer 2 ,eDielec2,Conduct2); 
HFieldPertbZCiy,ix)  =  FIFieldPertbTempCl,3) ; 

end 

end 


%  Find  absolute  Perturbed/Direct  (in  PPM)  for  FI  Field  for  each  of  the 
%  measurement  position  points  on  the  grid 
HFieldRatioPPM  =  1 . 0E6 . *absCHFieldPertbZ . /HFieldDirZ)  ; 

%  Plot  Perturbed  FI  Field  in  PPM  for  each  of  the  measurement  position 

%  points  on  the  grid 

figure 

pcolor(xx  ,yy ,  FIFieldRatioPPM) 
shading  interp 
xlabelC'x  -  meters') 
ylabelC'y  -  meters') 
colorbar 

titleStrLl  =  sprintf  ( '  Flz(perturb)/FlzCdi  rect)  [PPM]\n'); 
titleStrLZ  =  'A  0.5  m  Radius  Iron  Sphere  at  1.5  m  below  surface'; 
titleStr  =  [titleStrLl, titleStrLZ] ; 
titleCtitleStr)  ; 
axis  image 

^♦**************************************************************************** 

We  have  supplied  this  MATLAB  code,  LoopSensor Example. m,  along  with 
getMagDipMomCond . m  on  the  CD-ROM. 


The  results  of  this  calculation  are  shown  below  as  a  color  intensity  plot  of  fhe  ratio 

I  J_J  Perturb  ^  J-J  Direct  | 

for  fhe  grid  in  PPM;  fhe  resulfs  have  been  smoofhed  using  fhe  MATLAB  shading 
interp  command. 
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Menu  Interface  Routines 

This  section  describes  the  MATLAB  interface  that  allows  a  user  to  define  fhe  sources, 
receivers  and  buried  structures.  The  user  can  then  setup  appropriate  linkages  to  any 
computations  to  configure  an  easy  to  use  menu  driven  interface.  The  menu  interface 
software  has  been  developed  so  that  linkages  to  any  user  defined  MATLAB  code 
(including  fhe  Sommerfeld  routines  or  the  Hill /King  buried  wire  routines)  may  be 
easily  implemented  by  an  experienced  MATLAB  user. 

Objects 

1.  There  are  objects  called  "sources"  and  "receivers"  that  have  locations  that  vary;  i.e., 
they  can  move.  These  locations  are  called  "measurement  points"  or  "time  steps". 
However,  the  order  of  the  measurement  points  is  not  important;  they  are  not  time 
related.  There  may  be  any  number,  up  to  some  maximum,  of  measurement  points. 
There  must  be  at  least  one  measurement  point. 

There  may  be  any  number,  up  to  some  maximum,  of  fhe  objects.  For  example,  you 
may  have  1  source  and  3  receivers.  All  sources  and  receivers  have  the  same  number  of 
measurement  points.  There  must  be  at  least  one  source  and  one  receiver. 

2.  There  are  also  objects  known  as  "buried  structures"  that  have  locations  that  do  not 
vary;  i.e.,  they  cannot  move.  Again  there  may  be  any  number,  up  to  some  maximum, 
of  fhese  objects.  And,  there  may  be  no  (0)  buried  structures. 

3.  A  location  is  specified  in  Cartesian  coordinates  as  a  point,  (x,  y,  z),  in  meters 
referenced  from  an  arbitrary  origin,  (0,  0,  0).  A  measurement  point  is  also  specified  in 
meters,  (x,  y,  z),  from  fhe  same  arbifrary  origin. 

4.  All  sources  have  a  strengfh  fhat  is  given  in  units  appropriate  to  the  type  of  the  source 
There  are  four  fypes  of  sources:  (1)  loop,  (2)  electric  dipole,  (3)  magnetic  dipole  and  (4) 
straight  wire. 

a.  A  loop  source  is  defined  by  its  radius,  in  meters,  and  its  orientation,  the  plane  that  is 
perpendicular  to  a  specified  unit  vector  located  at  its  center  point.  The  unit  of  strengfh 
of  a  loop  source  is  "Amp".  The  location  of  fhe  center  of  fhe  loop  will  define  fhe  location 
of  fhe  loop.  The  loop  intended  here  has  finite  size  with  respect  to  other  dimensions  of 
fhe  problem,  e.g.  layer  depths,  distances  to  objects,  etc. 
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b.  An  electric  dipole  source  is  located  at  a  single  point  and  is  defined  by  its  orientation,  a 
specified  unif  vecfor  af  fhaf  poinf.  The  unif  of  sfrengfh  of  an  elecfric  dipole  source  is 
"Amp-m".  If  is  an  infinifesimal-lengfh  wire  fed  af  ifs  cenfer,  in  ferms  of  ifs  radiafed  fields. 

c.  A  magnetic  dipole  source  is  also  locafed  af  a  single  poinf  and  defined  by  ifs 
orienfafion,  a  specified  unif  vecfor  af  fhaf  poinf.  The  unif  of  sfrengfh  of  a  magnetic 
dipole  source  is  "Amp-m^".  If  is  an  infinifesimally  small  loop,  whose  axis  coincides  wifh 
fhe  perpendicular  fo  fhe  loop  cenfer,  wifh  respecf  fo  ifs  radiafed  fields. 

d.  A  sfraighf  wire  source  is  defined  by  ifs  lengfh  and  direcfion,  bofh  given  by  a 
specified  vecfor  wifh  cenfer  af  fhe  midpoinf  of  fhe  wire.  The  locafion  of  fhe  cenfer  of 
fhis  vecfor  will  define  fhe  locafion  of  fhe  wire.  The  unif  of  sfrengfh  of  a  wire  source  is 
"Volf".  If  fhe  wire  lengfh  is  shorf  (approaching  intinifesimal),  fhe  wire  is  usually  better 
or  more  convenienfly  modeled  by  fhe  elecfric  dipole.  The  feed  of  fhe  wire  is  assumed  fo 
be  af  fhe  cenfer  (fhis  may  be  changed  by  a  MATLAB  modification  fo  fhe  menu  inferface 
program). 

5.  There  are  four  fypes  of  receivers,  fhe  same  as  fhe  sources.  However,  receivers  do  nof 
have  an  associafed  sfrengfh.  All  ofher  properties  are  fhe  same  as  fhe  four 
corresponding  fypes  of  sources. 

6.  There  are  fwo  fypes  of  buried  sfrucfures:  (1)  ellipsoid  and  (2)  wire. 

a.  An  ellipsoid  buried  sfrucfure  is  defined  by  ifs  3  semi-axes,  ax,  ay  and  az;  ifs 
orienfafion,  which  may  be  specified  by  fwo  unif  vecfors  (one  corresponding  fo  ifs  local 
z-axis  and  one  corresponding  fo  ifs  local  x-axis)  or  fhree  rofafions  (angular  rofafions 
abouf  each  of  ifs  local  axes);  and  ifs  elecfrical  properties  which  are  ifs  relative 
permeabilify,  relafive  dielecfric  consfanf  or  permiffivify,  and  whefher  or  nof  if  is  a 
perfecf  conducfor.  The  cenfer  of  fhe  ellipsoid  will  define  ifs  locafion. 

b.  A  sfraighf  wire  buried  sfrucfure  is  defined  by  ifs  lengfh  and  direcfion,  bofh  given  by 
a  specified  vecfor  wifh  cenfer  af  fhe  midpoinf  of  fhe  wire.  The  cenfer  of  fhis  vecfor  will 
define  ifs  locafion. 

Global  Variables 

The  above  information  is  communicafed  fo  fhe  resf  of  fhe  program  via  fhe  following 
global  variables,  found  in  fhe  m-file  "SRLGlobals.m". 
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Sources  /  Receivers 

nMeasure  -  the  number  of  measurement  points. 
nSource  -  the  number  of  sources. 
nReceiver  -  fhe  number  of  receivers. 

xSource,  ySource,  zSource  -  arrays  (nMeasure,nSource). 

These  arrays  contain  the  measurement  locations  in  meters  of  all  sources. 

sourceType  -  array  (l,nSource). 

This  array  contains  the  type  of  each  source,  1  ==  loop,  2  ==  electric  dipole, 
3  ==  magnetic  dipole  and  4  ==  wire. 

sourceStrength  -  array  (l,nSource). 

This  array  contains  the  strength  of  each  source. 

sourceProperties  -  array  (nSource,4). 

This  array  confains  fhe  defining  properties  of  each  source.  There  is  one 
row  per  source  and  the  contents  of  each  row  varies  wifh  fhe  type  of 
source  as  follows: 

loop  -  columns  1-3  ==  fhe  orientation  unit  vector, 

column  4  ==  radius  of  fhe  loop  in  mefers; 
dipoles  -  columns  1-3  ==  fhe  orientation  unit  vector, 
column  4  ==  not  used; 

wire  -  columns  1-3  ==  the  orientation  and  length  vector, 

column  4  ==  not  used; 

xReceiver,  yReceiver,  zReceiver  -  arrays(nMeasure,nReceiver). 

These  arrays  contain  the  measurement  locations  in  meters  of  all  receivers. 

receiverType  -  array  (l,nReceiver). 

This  array  is  similar  to  the  sourceType  array. 

receiverProperties  -  array  (nReceiver,4). 

This  array  is  similar  to  the  sourceProperties  array. 
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Additional  source  /  receiver  variables 

The  following  variables  are  created  if  the  user  decides  to  have  measurement  locations 
computed  from  a  grid. 


xGrid,  yGrid,  zGrid  -  arrays  (4). 

These  arrays  contain  the  definition  of  fhe  computed  grid.  The  four 
elements  of  each  array  define  fhe  minimum  value  in  meters,  the 
maximum  value  in  meters,  the  number  of  grid  points  and  the  distance 
between  grid  points  in  meters  -  for  each  axis. 

xOSource,  yOSource,  zOSource  -  arrays(nSource). 

These  arrays  confain  fhe  location  of  each  source  relative  to  the  origin  of 
fhe  grid. 

xOReceiver,  yOReceiver,  zOReceiver  -  arrays(nReceiver). 

These  arrays  contain  the  location  of  each  receiver  relative  to  the  origin  of 
fhe  grid. 

Buried  Structures  (Ellipsoids) 

nElip  -  the  number  of  ellipsoids  defined. 

xElip,  yElip,  zElip  -  arrays  (l,nElip). 

These  arrays  give  fhe  locations,  in  meters,  of  fhe  centers  of  fhe  various 
buried  ellipsoids. 

ixElip,  iyElip,  izElip  -  arrays  (l,nElip). 
jxElip,  jyElip,  jzElip  -  arrays  (l,nElip). 
kxElip,  kyElip,  kzElip  -  arrays  (l,nElip). 

These  arrays  describe  fhree  unit  vectors  for  each  ellipsoid.  The  unit 
vectors  give  the  orientation  of  each  ellipsoid.  The  unit  vector  (jxElip, 
jyElip,  jzElip)  is  computed  by  the  program  given  the  other  two. 

dxElip,  dyElip,  dzElip  -  arrays  (l,nElip). 

These  arrays  give  the  lengths,  in  meters,  of  each  of  the  three  semi-axes  of 
each  ellipsoid. 
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pcElip  -  array  (l,nElip). 

This  array  provides  a  "perfect  conductor"  flag  for  each  ellipsoid.  (==  0 
means  the  ellipsoid  is  not  a  perfect  conductor). 

epsElip  -  complex  array  (l,nElip). 

This  array  provides  the  relative  complex  dielectric  constant  for  each 
ellipsoid. 

muElip  -  complex  array  (l,nElip). 

This  array  provides  the  relative  complex  permeability  of  each  ellipsoid. 

Buried  Structures  (Wires) 

nWire  -  the  number  of  wires  defined. 

xWire,  yWire,  zWire  -  arrays  (l,nWire). 

These  arrays  give  the  locations,  in  meters,  of  the  centers  of  the  various 
buried  wires. 

vxWire,  vyWire,  vzWire  -  arrays  (l,nWire). 

These  arrays  describe  a  vector  for  each  wire.  The  center  of  the  vector  is 
located  at  the  center  of  its  corresponding  wire.  And  the  direction  of  the 
vector  is  the  direction  of  the  wire.  And  the  length  of  the  vector  is  the 
length  of  the  wire  (meters). 

Private  Global  Variables 

The  source  /  receiver  and  buried  structures  information  is  organized  into  arrays  of 
structures  for  the  use  of  the  input  portion  of  the  program.  These  arrays  of  structures 
are  defined  within  a  "private"  global  area.  (See  the  m-file  "SRLPrivates.m".)  Here  are 
descriptions  of  these  private  global  variables. 

Sources  /  Receivers 

SRLmethod  -  the  method  of  input,  0  ==  by  computation,  1  ==  by  file. 

SRLpathname  -  the  path  and  name  of  the  input  file. 
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This  variable  is  used  if  the  input  is  by  file. 

SRLgridX,  SRLgridY,  SRLgridZ 

These  arrays  hold  the  same  information  as  the  variables  xGrid,  yGrid  and 
zGrid  described  above. 

SRLSn  -  the  number  of  sources. 

SRLSmax  -  the  maximum  allowed  number  of  sources. 

SRLSrc  -  An  array  of  structures  (a  cell  array). 

Each  structure  within  this  array  contains  the  relative  measurement 
location  and  defining  properties  of  one  source.  There  is  one 
structure  per  source.  The  cell  array  is  preallocated  with  the 
MATLAB  statement  "SRLSrc  =  cell(SRLSmax,l);". 

Here  are  the  elements  of  each  structure: 

.type  =  integer  -  the  type  of  source, 

1  ==  loop 

2  ==  electric  dipole 

3  ==  magnetic  dipole 

4  ==  wire 

.location  =  array  (1,3) 

The  location  (x,y,z  -  meters)  of  this  source  relative  to 
the  current  measurement  point 
.vector  =  array  (1,3)  the  orientation  (and/ or  length)  vector 
.radius  =  the  radius  of  a  loop  source. 

.strength  =  complex  -  the  source  strength; 

Specific  quantities  within  the  array  of  structures  are  referred  to 
with  notation  in  the  form  of  "SRLSrc{5}.location(2)",  which  refers  to 
the  y-coordinate  of  the  relative  location  of  the  5th  source. 

SRLRn  -  the  number  of  receivers. 

SRLRmax  -  the  maximum  allowed  number  of  receivers. 
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SRLRcv  -  this  is  a  cell  array  of  structures  that  define  fhe  locations  and 
of  fhe  receivers.  It  is  defined  similar  fo  fhe  above  cell  array, 
difference  is  that  the  value  of  fhe  "strength"  member  of  fhe 
wifhin  fhe  array  will  be  sef  to  zero. 

Buried  Structures 

BSDpathname  -  the  path  and  name  of  the  input  file. 

This  variable  is  used  if  fhe  input  is  by  file. 

BSDEn  -  fhe  number  of  ellipsoids. 

BSDEmax  -  fhe  maximum  allowed  number  of  ellipsoids. 

BSDEllip  -  cell  array  of  ellipsoid  buried  sfructures 

One  structure  per  ellipsoid  defined  as  follows: 

.location  =  array  (1,3),  the  location  of  fhe  cenfer  of  fhe  ellipsoid  in 
Carfesian  coordinates  (x,y,z  -  meters) 

.axes  =  array  (1,3),  the  lengths  of  fhe  semi-axis,  x,  y  and  z,  in 
meters 

.xunit=  array  (1,3),  unit  vector  describing  the  orientation  of  fhe 
ellipsoid's  x-axis  relative  to  its  local  origin 
.yunit=  array  (1,3),  unit  vector  describing  the  orientation  of  fhe 
ellipsoid's  y-axis  (this  vector  is  calculated) 

.zunit=  array  (1,3),  unit  vector  describing  the  orientation  of  fhe 
ellipsoid's  z-axis  relative  to  its  local  origin 
.angles  =  array  (1,3),  alternate  orientation  expressed  as  angles  of 
rotation  about  the  respective  axes  (degrees),  (the  x,y,z-unit 
vectors  are  computed  from  fhis), 

.orient  =  integer,  orientation  flag,  l==unit  vectors,  0==angles, 

.pc  =  integer,  perfect  conductor  flag,  l==yes,  0==no, 

.eps=  complex,  relative  dielectric  constant, 

.mu  =  complex,  relative  permeability  constant. 

BSDWn  -  the  number  of  wires. 


properties 
The  only 
structures 
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BSDWmax  -  the  maximum  allowed  number  of  wires. 


BSDWire  -  cell  array  of  wire  buried  structures 

One  structure  per  wire  defined  as  follows: 

.location  =  array(l,3),  the  location  of  fhe  midpoint  of  fhe  wire  in 
Carfesian  coordinafes  (x,y,z  -  meters) 

.vector  =  array(l,3),  a  vector  with  its  center  at  the  midpoint  of  fhe 
wire.  This  vecfor  describes  lengfh  (mefers)  and  direction. 


Operation 

To  input  source  and  receiver  locations  invoke  the  MATLAB  function  "SRLocations".  To 
input  the  buried  structures  definitions  invoke  the  MATLAB  function  "BSDefinition". 

Each  of  fhese  functions  will  present  a  dialog  box  that  allows  the  user  to  select  between 
input  by  computation  or  input  by  file. 

Input  via  File 

If  you  intend  to  select  the  input  by  file  mefhod  you  must  have  on  the  disk  a  file  fhat  is 
readable  via  fhe  MATLAB  "load  filename"  sfatement. 

If  your  file  inpuf  is  for  fhe  sources  and  receivers  then  your  input  file  must  have  all  of  fhe 
source  and  receivers  variables  described  in  fhe  sections  "Sources /Receivers"  and 
"Additional  source/ receiver  variables"  in  the  "Globals"  section  above. 

If  your  file  inpuf  is  for  fhe  buried  sfrucfures  fhen  your  input  file  must  have  all  of  fhe 
buried  structures'  variables  described  in  the  sections  "Buried  Structures  (Ellipsoids)"  and 
"Buried  Structures  (Wires)"  in  the  "Globals"  section  above. 

Source  /  Receiver  Input  via  Computation 

When  you  input  your  source/ receiver  definitions  via  computation  you  must  define  fhe 
grid  of  measurement  locations,  define  at  least  one  source  and  define  at  least  one 
receiver.  At  the  conclusion  of  your  input  the  definitions  will  be  converted  and  copied  to 
the  global  variables  described  in  the  "Globals"  section  above.  These  values  will  also  be 
written  to  the  file  "SRLSave.mat"  in  a  form  thaf  is  readable  via  fhe  MATLAB  "load" 
statement. 
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Each  time  the  "SRLocations"  function  is  invoked  it  will  look  for  fhe  "SRLSave.mat"  file 
and  read  it  and  present  those  values  to  you  when  you  select  input  by  computation.  If 
you  want  to  start  from  scratch  discard  the  file  "SRLSave.mat". 

The  dialog  boxes  presented  during  the  input  by  computation  are  self  explanatory.  Each 
box  represents  one  source  or  one  receiver.  You  navigate  from  source  fo  source  or 
receiver  to  receiver  via  the  buttons  within  the  box.  You  can  choose  to  add  or  delete 
items.  When  you  delete  an  item  it  is  moved  to  the  end  of  fhe  lisf  and  fhe  list  is 
shortened  by  one.  Thus,  when  you  subsequently  add  a  new  item,  the  one  you  deleted 
appears  at  the  end  of  fhe  list  for  user  information. 

Buried  Structures  Input  via  Computation 

When  you  input  your  buried  structures  definitions  via  computation  the  process  is 
similar  to  that  described  above.  At  the  conclusion  of  your  input  the  definitions  will  be 
converted  and  copied  to  the  global  variables  described  in  the  "Globals"  section  above. 
And  these  values  will  also  be  written  to  the  file  "BSDSave.mat"  in  a  form  thaf  is  readable 
via  fhe  MATLAB  "load"  statement. 

Each  time  the  "BSDefinition"  function  is  invoked  it  will  look  for  fhe  "BSDSave.mat"  file 
and  read  it  and  present  those  values  to  you  when  you  select  input  by  computation. 

And,  if  you  want  to  start  from  scratch  then  discard  the  file  "BSDSave.mat". 

The  dialog  boxes  presented  for  buried  sfructure  definitions  work  the  same  as  those  for 
source  /  receiver  definitions.  There  is  one  box  for  each  ifem.  You  may  have  zero  items, 
in  which  case  all  entries  in  the  box  are  disabled  except  the  "add"  button. 
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Hill/King  Buried  Wire  Routines 

For  a  number  of  applications  it  is  useful  to  model  long  insulated  conductor  in  a  lossy 
medium  with  arbitrary  termination  on  either  end.  We  have  written  MATLAB  routines 
that  use  David  A.  Hill's  paper  (expanding  on  previous  work  of  R.W.P.  King  -  "Magnetic 
Dipole  Excitation  of  an  Insulated  Conductor  of  Finite  Length",  IEEE  Trans.  Geosci. 
Remote  Sensing,  vol.  28,  pp.  289-294,1990)  to  do  these  calculations.  The  Hill /King  result 
allows  current  on  an  insulated  conductor  in  a  lossy  medium  with  arbitrary  termination 
on  either  end  to  be  obtained  from  a  defined,  imposed,  external  electric  field. 

The  Hill /King  model  assumes  that  the  wire  is  buried  in  a  single  infinite  medium;  we 
are  interested  in  multilayered  applications  -  so  there  is  a  issue  of  how  the  current  is 
affected  by  getting  near  a  medium  interface;  we  have  used  NEC-3  to  address  this  issue 
(discussed  in  detail  below).  It  is  our  opinion  that  the  single  infinite  medium 
approximation  for  the  calculation  of  the  Hill  /  King  current  will  work  well  for  multilayer 
applications.  Included  on  the  CD-ROM  are  examples  using  this  code;  the  examples  are 
taken  from  Hill's  paper  reproducing  Figures  6-10  in  the  paper. 

To  determine  the  effect  of  a  medium  interface  on  a  current  distribution  we  used  NEC-3 
to  look  at  the  current  distribution  as  a  voltage  excited  horizontal  buried  wire  moves 
closer  to  a  medium  interface.  We  used  3  kHz  and  average  earth  (s=10,  a=0.01  mhos/ m). 
We  looked  at  two  lengths  of  wire;  one  a  little  under  2.5  A.n,(,dmm  in  length  (1000  m)  and 
one  a  little  under  0.5  k^edium  in  length  (200  m).  The  wire  was  fed  at  the  center  with  1  Volt. 
The  depths  used  were  2000  m,  200  m,  20  m,  2  m,  and  0.2  m.  The  results  of  these 
calculations  show  that  there  is  a  small  change  in  the  magnitude  of  the  current  and  a 
slight  increase  in  the  phase  as  the  wire  approaches  the  interface. 


ELF/VLF  SBIR  Phase  II  Final  Report 


CODAR  Ocean  Sensors,  Ftd. 


Page  21  of  42 


Current  Magnitude  for  Center  Fed,  Horizontal,  Buried  Wire  at  Various  Depths  Using  NEC-3  at  3  kHz 
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Current  Phase  for  Center  Fed,  Horizontal,  Buried  Wire  at  Various  Depths  Using  NEC-3  at  3  kHz 
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Current  Magnitude  for  Center  Fed,  Horizontal,  Buried  Wire  at  Various  Depths  Using  NEC-3  at  3  kHz 
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Position  along  wire  [m] 


Current  Phase  for  Center  Fed,  Horizontal,  Buried  Wire  at  Various  Depths  Using  NEC-3  at  3  kHz 


o 


00  CO  (M  CM  1- 

[Bea]  iLuejjnQ  jo  esegd 


ELF/VLF  SBIR  Phase  II  Final  Report 


CODAR  Ocean  Sensors,  Ftd. 


Page  25  of  42 


Position  along  wire  [m] 


Appendix  A  -  Sommerfeld  Integral  Definitions 
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VED  in  Region  1 

Located  at  z  =  0  -  Interface  1-2  at  z  =  -d^  -  Interface  2-3  af  z  =  -d^ 


Fields  in  Region  1:  (z  ^  -d, ) 
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VMD  in  Region  1 

Located  at  z  =  0  -  Interface  1-2  at  z  =  -d^  -  Interface  2-3  af  z  =  -d^ 
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VED  in  Region  2 

Located  at  z  =  -  4  (-^2  ^  -^o  <  “^1)  "  Interface  1-2  at  z  =  -d^ 


Fields  in  Region  1:  (z  ^  -d^) 
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Fields  in  Region  2  above  the  source:  (-d,  >  z  ^  -Zo) 
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VMD  in  Region  2 

Located  at  z  =  -  4  (-^2  ^  -^o  <  “^1)  "  Interface  1-2  at  z  =  -d^ 
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HED  in  Region  1 

Located  at  z  =  0 

Fields  in  Region  1:  (z  >  -d, ) 


Interface  1-2  at  z  =  -dj 


Interface  2-3  af  z  =  -d. 
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Fields  in  Region  2:  (-d,  >  z  >  -d^ ) 
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HMD  in  Region  1 

Located  at  z  =  0 

Fields  in  Region  1:  (z  >  -d, ) 
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Interface  1-2  at  z  =  -dj 


Interface  2-3  af  z=  -d. 
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Fields  in  Region  2:  (-d,  >  z  >  -d^ ) 
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jyTM 
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HED  in  Region  2 

Located  at  z  =  -  4  -  Interface  1-2  at  z  =  - 

Fields  in  Region  1:  (z  ^  -d^) 
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Fields  in  Region  2  below  the  source:  (-4  >  z  >  -^2 ) 

n  ”  ,  .  , 


pH£D[2]  _  ■  II 

^2Z  ~  ^  A 

Ajtcoe. 


00 

cos4>^dkpklJ,[kpp)d. 

0 

71  iP' 

-  i—im4,pk^-^j{k,p]c 


EMl^-ik2zZ  _!_  d^iz  (z+2d2)j^TM 


TeI  ^-ikzz  z  _!_  d'^iz  {z+  2d2  TE 


T^HEDil^ 


—msUdkJcz-kf  J<,(kiJ))- 

2  0  I,  kpP  j 


^EM  _^~^^2zZ  ^^k2z{z  +  2d2) j^TM 


-£ 


II 


-cos 


(j)J  dk^ 


0  K 


^27 


r7HEE(2^ 

Ejzh  ^  ^  = 

^  Ajiooe^p 


Anooe^p 

71  00 

0 


+kl 


II 


Ajiooe. 


0  ^2z\  ^pP 
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4jrp 

■£™'^/‘'^a[-^.(*pp)-  (-«*'= +e'‘“'”'‘’4^ 

0  \  P^  ) 

H,"r'  - 

n 

-  ^cos  0J<iA:/,(A:^p)cr(-e-“-  +  "  “■  >RS] 
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HMD  in  Region  2 

Located  at  z=  -  4  -  Interface  1-2  at  z=  -dj  -  Interface  2-3  af  z=  -d^ 


Fields  in  Region  1:  (z  >  -d^) 

An  £1  J  ^^2^ 

-  -  ■^^^cos  Lpy*"'""'  'e-‘ 

4yr^i  Q 


21  ^^-^2 


±21  iyi2 


Anp 

/w/TjM  £2 

4jr  £; 


jpTEii^TE 

±21  iyi2 


■Sin 


(Pfdk^^k^ 

0  '^2z 


Kp 


^i^\z  (^+'^1  ) ^-ikizdirpTM 


^hmd[2]  ^  ^  COS  (pjdkpk pi  Jo(kpp)  - 


■^l{^pp)\ 


kpP  j 


^ikiz(z+di)^-ik2^dy  rYTEj^TE 


icopi^lA  82 

4jtp  8^ 


7 

coscpjdkp  !h^J^i^kpp)e‘'''^^^"‘‘'  fj 


Ml 


-0  K 


HMD[2]  ^ 


‘IP 


+k 


Anp^ 

2  ip2l^  ^2 


J,(Kp) 
k,P  ! 


Anp^p  £ 


00  j 

(l>Jdkpk^^kp^Jo(kpp)- 

GO 

COS  72 


21  "■‘2 


^/lij(z  +  c(i  )^-ik2Pi'pTE  j^TE 


M 

21  ^"2 


1  0 


^HMD[2]  ^ 

‘\np^p  •' 


-ik2pt  kfzTE  llyfTE 


721  Ml 


” 


-0  ^2^ 


J,(k,p)' 
kfP  , 


^ikidz+d^  )^-ik2zdi'pTM j^TM 


Fields  in  Region  2  above  fhe  source:  (-d,  >  z  >  -Zq) 

0  ^2z 

T  A 

^HMD[2]  ^  _  cos(l)JdkpklAkpp)M2^i^e‘’'^^^  +  e 


ikz^z  ^  ik2  2(z+2d^  ) j^TM 


-ik2z(zA2d,)  TE 
^21 


pHMD[21  i(t)P2lA 

~  Anp 


sm(l)JdkpM^^J^[kpp)l. 


ik2zZ  ^  „-ik2z{z  +  2dd j^TE 


21 


+  sin  (t)fdkpM™kp 

0 


4jr 


'^o(^pP) 


\ 


Kp 


^ikzzZ  _  ^-‘k2z(z+2d d j^M 
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^hmd[2]  ^  COS (pJdkpM^^kp  Jo(kpp)- 


4kp) 


{z+2di)j^TE 


H“‘^^l--—cos^Jdk^t,Mn^  J„(k^p)- 


J,{k,p) 


{z+2di)j^TE 


-kl  ^cos^Jrf*^  +  e-"-!”"'  Ij!” 


'T  \  ^ 


-/^2  J  (  ^  ■'■^  ^1  )  D 


g'^Z  _|_  g-*^2!(z+2«'i  )^7M 


Fields  in  Region  2  below  the  source:  (-4  >  z  >  -dj ) 

,  ,  CO  72 


■iljjZ  .  Jl^2z{z  +  2d2)  nTM 


HzT^'--^^o^<PfdkX4KpK{‘ 
E;r"'  ■  ‘■^^^iffpd^pKi'^ 

+  “°^^'^sin0jrfA:/p  2„(j:,p)- 


“'^Z  ^  (z+2d2)j^TE 


^pP 

2i(v) 


/_g  -  '*2zZ  _j_  g*^2z  (z+  2^2  ) 


/a)p2^d 


cos0jd/t/i(/tpp)iV2™(-c-'^^^  +  c' 
0 


, ^^2  zi.^ ’’’^'^ 2  ) 

^3 


^™d[2]  ^  ^  j  _ 


TeI  -ik2zZ  .  ik2z{z+2d2)  nTE 


-e-^zzz  +  ^ 


00 

-k^,  ^cos^fdk^  ±J^(k^p)Nr{e-‘^^^^  + 


!i:2  j(z+2  4  )  pTM 
^3 
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^ik2-  {z+2d2)  pTE 
^23 


^k^l^sm(l)Cdk^-^  Jo(kpp)- 

4jr  K2^ 


J,(Kp]' 

kpP 


J\J™le-ik2^Z  +  JI^TM 
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Appendix  B  -  Contents  of  CD-ROM 

On  the  CD-ROM,  named  ELF_VLF  delivered  with  this  report  are  three 
folders/  subdirectories: 

Hill_King,  Menujnterface,  Sommerfeld_Routines 

In  addition  to  these  folders/ subdirectories  is  a  PDF  file  version  of  this  Report 

ELF_Report.PDF” 

Below  are  listed  the  files  contained  in  each  of  the  folders/ subdirectories 

Hill_King: 

HillFiglO.m 

HillFige.m 

HillFigZ.m 

HillFigB.m 

HillFigg.m 

KingsGreenFtn.m 

Menujnterface: 

BSDefinition.m 

BSDefinitionBox.m 

BSDefinitionBox.mat 

BSDEIlipsoid.m 

BSDEIlipsoidBox.m 

BSDEIlipsoidBox.mat 

BSDWireStruct.m 

BSDWireStructBox.m 

BSDWireStructBox.mat 

ellipplot.m 

SRLGIobals.m 

SRLGrid.m 

SRLGridBox.m 

SRLGridBox.mat 

SRLObject.m 

SRLObjectBox.m 

SRLObjectBox.mat 

SRLocations.m 

SRLocationsBox.m 

SRLocationsBox.mat 

SRLPrivates.m 

Sommerfeld_Routines: 

EFIEDRegion1_1  .m 
EFIEDRegion1_2.m 
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EHEDRegion2_1  .m 

EHEDRegion2_2.m 

EHMDRegion1_1  .m 

EHMDRegion1_2.m 

EHMDRegion2_1  .m 

EHMDRegion2_2.m 

EVEDRegion1_1.m 

EVEDRegion1_2.m 

EVEDRegion2_1  .m 

EVEDRegion2_2.m 

EVMDRegion1_1  .m 

EVM  D  Reg  ion  1  _2 .  m 

EVMDRegion2_1  .m 

EVM  D  Reg  ion  2_2 .  m 

ExampleHED.m 

fastint.m 

fastIntSkIz.m 

getElecDipoleField.m 

getkrinterp.m 

getMagDipMomCond.m 

getMagDipoleField.m 

FIFIEDRegion1_1  .m 

FIFIEDRegion1_2.m 

FIFIEDRegion2_1  .m 

FIFIEDRegion2_2.m 

FIFIMDRegion1_1  .m 

FIFIMDRegion1_2.m 

HHMDRegion2_1  .m 

HHMDRegion2_2.m 

FIVEDRegion1_1.m 

FIVEDRegion1_2.m 

FIVEDRegion2_1  .m 

FIVEDRegion2_2.m 

FIVMDRegion1_1  .m 

HVMDRegion1_2.m 

FIVMDRegion2_1  .m 

HVMDRegion2_2.m 

interpICOS.m 

LoopSensorExample.m 

psiN.m 

SommerfeldGlobals.m 

Sommerfeldinit.m 
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